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Abstract 



The fully connected -fC-user interference channel is studied in a multipath environ- 



ment with bandwidth W. We show that when each link consists of D physical paths, 
the total spectral efficiency can grow linearly with K. This result holds not merely in 
the limit of large transmit power P, but for any fixed P, and is therefore a stronger 
characterization than degrees of freedom. It is achieved via a form of interference 
alignment in the time domain. A caveat of this result is that W must grow with K, 
a phenomenon we refer to as bandwidth scaling. Our insight comes from examining 
^ . channels with single path links {D = 1), which we refer to as line-of-sight (LOS) 

links. For such channels we build a time-indexed interference graph and associate the 
communication problem with finding its maximal independent set. This graph has 
^ ■ a stationarity property that we exploit to solve the problem efficiently via dynamic 

•/^ i programming. Additionally, the interference graph enables us to demonstrate the ne- 

I cessity of bandwidth scaling for any scheme operating over LOS interference channels. 

' Bandwidth scaling is then shown to also be a necessary ingredient for interference 

Q>^ . alignment in the K-usei interference channel. 

o 

oo 

^ : 1 Introduction 

> ■ 

^ ■ The problem of communicating efficiently in wireless adhoc networks has received much 
. attention of late, the focus being on how to deal with interference in a shared medium. 
Traditional approaches based on orthogonalizing users (eg. TDMA or FDMA) or reusing 
the spectrum (eg. CDMA, certain modes of 802.11) suffer from poor spectral efficiency. In 
particular as the number of users in the system grows, the spectral efficiency of each link 
vanishes. More recent approaches include using multi-hop P, distributed MIMO [2], [3], 
and interference alignment [1], [5], [6]. This work belongs to the latter category. Interference 
alignment is a technique that uses appropriate precoding to compact interfering signals into 
small dimensional subspaces at each receiver. At the same time, the subspace occupied by 
the data remains linearly independent of the interference. It was first applied to a multiple 
MIMO base station problem in [1] and shown to be capable of achieving multiplexing 
gains distinctly greater than those achievable using conventional signalling schemes. The 
technique was then extended in [6] to show that there were exactly 4M/3 degrees of freedom 
in the MIMO X channel (two MIMO transmitters each desiring to send data to two MIMO 
receivers) with M > 1 antennas at each transceiver. Following this, a more sophisticated 
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interference alignment technique was developed in [5] for the i^-user interference channel 
with an infinite number of independently faded sub-channels, and used to demonstrate that 
contrary to conventional wisdom, the total degrees of freedom of the channel is K/2. 

Whilst the last approach demonstrates the potential benefits that interference alignment 
techniques can provide, it suffers from a number of limitations. Perhaps the foremost is 
that whilst a degrees of freedom characterization is useful in the high SNR limit, it may 
not be meaningful at moderate SNR's. This stems from the fact that degrees of freedom 
characterizes the asymptotic slope of the spectral efficiency curve and not its actual value. 
In particular, it is unclear whether at any fixed SNR the total spectral efficiency of the 
system is increasing in proportion to K, or increasing at all. The point here is that [5] 
does not contain a scaling law result the likes of [1] and [2], that is, it does not tell us what 
happens to the system capacity as more users enter the fray. We address this question by 
constructing a communication scheme that achieves a scaling of system capacity arbitrarily 
close to linear. 

There is another key limitation. It is natural to interpret the parallel channels used 
in the interference alignment scheme of [5] as corresponding to sub-channels in the fre- 
quency domain. This is due to various difficulties associated with realizing the scheme over 
independently faded parallel channels in time, most notably the very rapid and accurate 
channel measurement that must take place, and the substantial delay incurred. But in 
order for a large number of frequency channels to undergo independent fading, significant 
scattering/multipath is required. 

In this work we examine the i^'-user interference channel with limited multipath. We 
start by assuming each of the links consists only of a single physical path with complex 
gain hij and delay Tij seconds. This model is a good representation for a line-of-sight 
(LOS) channel. Following this we generalize to the case where each link consists of D 
physical paths. We illustrate a simple and elegant representation of interference alignment 
in the time domain, in terms of aligning symbols on a time-indexed interference graph. 
This interference graph proves to be an extremely useful tool for both conceptualizing 
and solving various problems relating to interference in LOS channels. We identify the 
problem of communicating on the LOS interference channel with the problem of finding 
a maximal independent set in the interference graph, and show how this problem can be 
solved efficiently using dynamic programming principles. The simplicity of this approach 
makes it quite versatile and potentially capable of being extended to tackle a variety of 
related problems. 

For the remainder of this work power spectral density (PSD) is used in place of signal- 
to-noise ratio, as we will later wish to compare schemes that use different bandwidths. 
Depending on the link delays, it may be possible to achieve a spectral efficiency as high 
as I log2(l + \hii\'^PSD / Nq) bps/Hz for each user i, regardless of the number of interferers. 
This is exactly half the spectral efficiency achievable in the absence of interference. We 
characterize the precise channel conditions for which this is possible and show that they 
occur in at least 1/16 of all scenarios for the 3-user interference channel. 

This of course, says little about the typical gains one can expect by aligning interference 
in the time domain. To address this question, we treat the link delays as independent 
and uniformly distributed random variables in [0,Td), where is the delay-spread of the 
channel. Previous work such as [1], [5], [6] has focused on degrees of freedom as a metric 
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for performance, which is a measure of the scahng of spectral efficiency with PSD when 
the power spectral density is large. We focus on the scaling of spectral efficiency with K 
when the number of users is large. The main result of this work is the construction of a 
communication scheme that achieves a spectral efficiency arbitrarily close to 0(1) log2(l + 
\hii\^PSD / Nq) for each user i a,s K oo. Thus our result characterizes the best possible 
scaling of spectral efficiency with K for any fixed PSD, as compared with [Ij, [5], [6], where 
the best possible scaling of spectral efficiency with PSD, for any fixed is characterized. 
In this sense, our characterization has a similar flavor to characterizations of scaling laws 
for wireless adhoc networks such as [1], however, our scheme requires no cooperation 
between users. 

A caveat of our result is that the bandwidth must scale sufficiently with K. Interpreting 
the parallel channel model of [5] in the frequency domain, one sees that the bandwidth 
there must also scale with K. Interestingly enough, the bandwidth scaling required for our 
scheme is essentially the same as that required in [5], namely 0{K^^ ). However, whereas 
the scheme in [5] requires coding over blocks of length 0{K^^ ), which creates significant 
encoding and decoding complexity issues as well as substantial delay, our scheme requires 
no block coding, and consequently does not incur any delay or suffer from complexity issues. 

Essentially, in order to align interference into a small dimensional subspace and keep 
it linearly independent of the data subspace, a high degree of resolvability of the received 
signals is required. We establish this concretely for the LOS channel in the context of 
the interference graph, showing that if the bandwidth scales sub-linearly with K, then the 
total spectral efficiency of the system (the sum of all users spectral efficiencies) will scale 
sublinearly, and hence almost all of the users will witness vanishing spectral efficiencies as 
K increases. 

This suggests that the greatest performance gains can be reaped in systems with large 
delay spreads. Perhaps the best example of such a system is the backbone of a mesh net- 
work, used to wirelessly connect rural areas. Such systems are well approximated by a 
LOS model, have large delay spreads, and are relatively static, making channel measure- 
ment simpler and more accurate. 

The structuring of the rest of the paper is as follows. In section [2] we describe the 
model of the i^-user LOS interference channel. Section [3] provides a summary of the main 
result concerning the achievability of non-vanishing spectral efficiencies as the number of 
users grows. The time-indexed interference graph is introduced in section HJ In the same 
section we present an algorithm for optimizing the spectral efficiency efficiently via dynamic 
programming. We also address the questions of bandwidth scaling, and of characterizing 
the class of channels for which the spectral efficiency can reach its maximum value. In 
section [5] we present our construction that establishes the main result. Following this, in 
section [HI we establish the relationship between time and frequency domain interference 
alignment techniques. Section [7] contains further discussion, extensions and open problems. 

2 Model 

We consider the i^'-user interference channel in which there are K transmitters and K 
receivers. Transmitter i wishes to send data to receiver i but its transmission constitutes 
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interference at all other receivers. We often refer to each tx-rx pair as a user. There are 
thus K"^ links in total, K direct links and K{K — 1) cross-links. Each link consists of a 
single physical path. Denote the gain and delay (in seconds) of the link between transmitter 
i and receiver i by hij G C and Tij G [0, cxd), respectively. We assume the hij and Tij are 
fixed for the duration of communication. Denote the signal transmitter j sends by Xj{t). 
Then the baseband signal at the zth receiver is 

K 

Viii) = hijXj{t - Tij) + Zi{t) 
i=i 

where the Zi{t) are i.i.d. white noise processes with power spectral density A'^q Watts/Hz. 
Denote the carrier frequency fc and the bandwidth that the signals Xj (t) are constrained to 
lie in by W Hz. Assuming the use of ideal sine pulses, the passband signal after sampling 
is given by 

K oo 

yi[m] = ^2 /iij-e"^"-^"'^'-' ^ sinc(/ - TijW)xj[m - I] + Zi[m] 
i=i z=o 

where the Zi[m] are i.i.d. CM{Q,NqW). We use the following conventional approximation 
for the sine pulse, 

sinc(t) ^ < t < 1/2; 

\ 0, otherwise, 

see page 27 of [7]. Let kj denote the integer round-off of the real number TijW. This leads 
to 

K 

yi[m] = ^ hije''^'^^''''''Xj[m - kj] + Zi[m]. 
i=i 

In a wireless model one typically makes an assumption about the statistics of the chan- 
nel. For example, in a channel with ISI, the tap coefficients are often modeled as i.i.d. 
Rayleigh random variables. In models where there is a dominant path, Rician random 
variables are used instead. Likewise in this work we make a statistical assumption on the 
channel, but limit this statistical assumption only to the link delays, taking the Tij to be 
i.i.d. uniform in [0,Trf), where denotes the delay-spread of the channel in seconds. This 
means that if we define L to be one plus the integer round off of T^W , the Uj are i.i.d. 
uniform in{0,...,L — 1}. No assumption is made on the link gains hij, other than they 
all being non-zero. 

We refer to this as the i^-user line-of-sight (LOS) interference channel, as this is the 
most common scenario giving rise to such a model. 

There is a straightforward extension of this model to the case where each link consists 
of D physical paths, such that the ith received signal is 

K D 

Vii^) = hij^dXjit - Tij^d) + Zi{t). 

i=i d=l 

Here hij^d G C and Tij^d G [0, oo) are the complex gain and delay of the dth. physical path 
between transmitter j and receiver i. This leads to the following generalization of the 
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passband model after sampling 

K D 

yi[m] = X] %,de"^^^=^»^-^Xj[m - kj^d] + Zi[m]. 

j=l d=l 

The natural extension of our statistical assumption for the LOS channel is to treat the 
delays Uj^d as i.i.d. uniform in {0, . . . , L — 1}. In doing so we are assuming independent 
delays not just across physical paths of different links, but also across the physical paths 
corresponding to the same link. 

We refer to this as the K-user D-path interference channel. 



3 Preview of Main Result 

Theorem 3.1. For any e > 0, there exists a communication scheme on the K-user LOS 
interference channel such that ifW > {2K{K — , the expected spectral efficiency 

of user i tends to 

1 / PSD\ 

as K ^ GO. 

Section O is devoted to proving this result. Here the expectation of spectral efficiency 
is taken over the random direct delays la. Roughly speaking this result says that as K 
scales, it is possible for each user to communicate at a spectral efficiency arbitrarily close 
to 0(1) log2(l + ^), so long as the bandwidth scales as fast as 0{{2K{K -l))^^^-^^). In 
other words, communication at spectral efficiencies that vanish arbitrarily slowly with K 
is possible if the bandwidth scales sufficiently. In this result, and throughout this work, we 
assume that all cross delays lij for i ^ j are known at all transmitters and receivers. This 
result is of a similar nature to the scaling laws of [2J , in that we show the growth of system 
capacity with the number of users is arbitrarily close to linear. Unfortunately, as is the 
case in |5], the required bandwidth scaling is great. A result presented later on (theorem 
14. 2p addresses the question of whether bandwidth scaling is necessary. 

For the case where each link consists of D physical paths we have the following gener- 
alization. 

Theorem 3.2. For any e > 0, there exists a communication scheme on the K-user D-path 
interference channel such that if W > {2DK{DK — i^^DK{DK-i)+e^ expected spectral 
efficiency of user i tends to 

1 , / ,, ,2PSD\ 

iDKiDK-l)y + ..S""^} l^-^l 1% ) 

as K ^ oo. 

As is evident from the statement of the theorem, there is a tradeoff between the number 
of physical paths per link and the bandwidth scaling required. Again, the question of 
whether this bandwidth scahng is necessary is addressed in a result presented later on. 
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t = t=l t = 2 t = 3 t = 4 t = 5 t = 6 t = 7 



Figure 1: The interference graph associated with the LOS channel with direct-delays In = 
0, 122 = 0, 133 = and cross-delays ki = 3, /31 = 1, lu = 1, ^32 = 4, /13 = 3, 123 = 0. 

4 The Interference Graph 

The key insight leading to theorem 13.11 comes from formulating the communication problem 
in a graph theoretical setting. We start with an example. Consider a 3-user LOS interfer- 
ence channel where the direct links are all zero, i.e. In = 0, I22 = and ^33 = 0, and the 
cross-links are say, I21 = 3, /31 = 1, I12 = 1, I32 = 4, /13 = 3, I23 = 0. Choose a length T for 
the communication block. Create a directed graph ^3,t(^i2, ^13, ^21, ^23, ^31, ^32) = as 
follows. Let the vertex set be 

V = MO), v,{T - 1)} U {V2{0), . . . , V2{T - 1)} U {V3{0), . . . , V3{T - 1)}. 

The vertex Vi{t) represents the tth time slot for the ith transmitter. Form the edge set 

5 as follows. Add a directed edge 621(0) starting from vertex fi(0) and ending at vertex 
'^2(^21) = 1^2(3). This represents the fact that, owing to a delay of 3 time slots, a transmission 
during time slot by transmitter 1 arrives at receiver 2 during time slot 3. Also add a 
directed edge 631(0) starting from fi(0) and ending at f3(/3i) = ^3(1). This represents the 
fact that, owing to a delay of 1 time slot, a transmission during time slot by transmitter 
1 arrives at receiver 3 during time slot 1. Likewise add directed edges 612(0) and 632(0) 
from vertex ^2(0) to fi(/i2) = "^1(1) and f3(/32) = "^3(4), respectively, and directed edges 613 
and 623 from vertex ^3(0) to fi(/i3) = fi(3) and f2(/23) = "^^2(0), respectively. This set of 
six edges encapsulates all of the interference generated by transmissions during time slot 0. 
As the channel is time-invariant the same interference structure applies for later time slots. 
Thus for each t = l,2, ...,T — 1 add a directed edge from Vi(t) to V2(t + /21), provided 
t + hi < T, a directed edge from vi(t) to vsit + I31), provided t + /31 < T, a directed edge 
from V2(t) to vi{t + I12), provided t + /12 < T, etc... See figure [T]for an illustration. 

In this example all direct delays were zero, whereas in general this is not the case. 
However as each transmitter j can merely offset it's transmitted sequence Xj[m] by —Ijj, 
we can effectively assume without loss of generality that Ijj = 0. More concretely we define 
the normalized cross- delays 

hj ''13 ■ 

Note that I'^j G {—L + 1,...,L — 1}, that is, it is possible for Z-^- to be negative. At this 
point one may wonder why the interference graph need be directed, since the feasibility of 
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etc... 



t = t=l t = 2 t = 3 t = 4 t = 5 t = 6 t = 7 

Figure 2: A feasible transmit pattern corresponds to an independent set. The vertices of a 
maximal independent set are shaded. For a symmetric channel, the maximal independent 
set maximizes total spectral efficiency. 

a given transmit pattern is independent of edge direction. The answer is it need not be, 
but we define it as such to aid in conceptualizing the problem. 
In general we have 

Definition 4.1. The time-indexed interference graph (or simply interference graph for 
short) of length T associated with the K-user LOS interference channel with normalized 
cross-delays {l'ij}i^j, is the directed graph GK,T{,{l'ij}i^j) — (V,^) where 

K 

V^[j{v^{l),...,v^{T)} 

i=i 

K K 

^-U U {e.,(l),...,e.,(T -/;.)}, 

with edge eij{t) stemming from vertex Vj(t) and ending at vertex Vi{t + 

This graph has KT vertices and approximately K{K — 1)T edges. As a transmission 
during time slot t is interfered with by one time slot at each other user, and itself interferes 
with one time slot at each other user, each vertex has both in-degree K — 1 and out-degree 
K — 1. To reduce the notational burden, we often refer to the graph QK,T{{l'ij}i^j) simply 
as Q. where the parameters of the interference graph are implicit. 

One can similarly define a time-indexed interference graph for the D-path interference 
channel, but in the interest of brevity and clarity we do not discuss it here. 

A transmit pattern is a subset of time slots during which data symbols are sent, one data 
symbol being sent per time slot. A transmit pattern is called feasible if each data symbol 
is received during a time slot that contains no interference from other transmissions. Thus 
a feasible transmit pattern corresponds to an independent set on the interference graph. 
Occasionally we will drop the adjective "feasible" when it is clear that the transmit pattern 
in question is such. As each data symbol arriving at receiver i during a time slot containing 
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no interfering symbol, is capable of conveying log2(l + bps/Hz, user z's spectral 

efficiency will be 



where iVj is the number of vertices in {fi(l), . • • ,Vi{T)} that are in the independent set. 
Let us assume for the meantime that ha = 1 for all i. Then the total spectral efficiency is 
directly proportional to the size of the independent set. Thus the problem of designing a 
communication scheme to maximize total spectral efficiency reduces to finding the maximal 
independent set of the interference graph. Denote the size of the maximal independent set 
of a graph Q (called the independence number) by q;(^). Then the maximum total spectral 
efficiency for a graph Q is 



For the preceding example, the maximal independent set is illustrated in figureO When- 
ever an independent set contains two vertices that possess a mutual neighbor, the interfer- 
ence generated by these two transmissions aligns at the mutual neighboring vertex. This 
is interference alignment in the time domain. For the example in figure [2l the neighbors of 
each unshaded vertex are all shaded, that is, data transmissions occur at all neighbors of an 
unshaded vertex. This is not always the case. Suppose we use a TDMA based communica- 
tion scheme where user 1 transmits on consecutive time slots for a long period whilst users 
2 and 3 remain silent. Then after a small guard interval user 2 transmits on consecutive 
time slots whilst users 1 and 3 remain silent. Finally user 3 transmits, and then back to 
user 1 and so on. In this round robin scheme, each unshaded vertex is connected to only a 
single shaded one and no interference alignment occurs. 

The problem of finding the maximal independent set is a well-known NP-hard problem, 
meaning that for an arbitrary graph, there is no known algorithm capable of solving the 
problem in time sub-exponential in the number of vertices. Knowing this it may appear 
that finding an optimal transmit pattern requires a computation time that is exponential in 
the block length T, however the interference graph is not an arbitrary graph. In particular 
it is stationary in the sense that, ignoring boundary effects, the structure of the graph is 
invariant to time shifts. In the next section we present an algorithm that exploits this 
property to find the maximum independent set in linear time. More generally, when the 
link gains ha are arbitrary the algorithm solves the problem of finding an independent set 
that maximizes spectral efficiency. We refer to this set as the optimal independent set. 

4.1 Finding the maximal independent set efficiently 

In this section we concentrate on the LOS channel, but the ideas can be extended to the D- 
path channel. Given an interference graph Q we now illustrate how dynamic programming 
principles can be employed to compute the maximal independent set efficiently. Let each 
vertex Vj{t) G {0, 1} with Vj(t) = 1 if Vj(t) is included in the transmit pattern, that is, if a 
data symbol is transmitted by user j during time slot t, and Vj(t) = otherwise. There is 
a slight abuse of notation here as we have used Vj{t) to represent both an element of the 
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t = t=l t = 2 t = 3 t = 4 t = 5 t = 6 t = 7 



Figure 3: The LOS interference channel of figure [H The state at time t = 4 is a function 
of the shaded vertices, each taking on one of two values. 



vertex set V and an indicator function for whether or not a data symbol is transmitted by 
transmitter j during time slot t. 

Concretely stated, the optimization problem we will solve is 

K /T-l \ 

min ^ii^)^ ] (3) 



s.t. Vi{s) + vj(t) < i,w(v^(s),vj{t)) ee 



.t=o 



where = log2(l + PSD /Nq). This is the problem of finding the optimal independent 
set. This cost is just the sum of the spectral efficiencies of the users weighted by the number 
of data symbols they send. We use a negative sign so as to justify the description of this 
metric as a cost, i.e. something we are trying to minimize. In the event that all direct 
gains are equal, rj is independent of j and the problem reduces to finding the maximal 
independent set of Q. 

To solve this problem efficiently, we start by defining 

I* = maxjmax/'., max— /'.}. 

This is the length of the longest edge that connects a vertex at a time t, to a vertex belonging 
to user i for time < t. When all the l[j are positive, /* simply represents the longest edge 
stemming from user j. In the example of figure [H = 3, = 4 and /g = 3. When some of 
the l[j are negative, /* represents the longest edge connecting user j to another user in the 
forward time direction. Thus the total amount of memory in the system is maxj /*. If this 
seems somewhat contrived, recall that although the interference graph is a directed graph, 
it need not be defined as such, as the effect of vertex Vi{t) causing interference at vertex 
Vj{t') is identical to the effect of vertex Vj{t') causing interference at vertex Viit). What 
matters for the dynamic programming formulation in this section, is not whether Vi{t) is 
causing interference with Vj{t') or vice versa, but whether t > t' , t = t' or t < t'. In our 
algorithm we move through vertices in order of increasing time t. The state of the system 
at time t is defined by those vertices at times < t that are connected to vertices at times 
> t. 
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State at t = 3 




t = t=l t = 2 t = 3 t = 4 t = 5 t = 6 t = 7 



Figure 4: A LOS interference channel with normahzed cross delays /21 = 0, /31 = 1, l'i2 = 2, 
Z32 = 0, /'i3 = 1 and = —2. The state at time t = 2 is a function of the shaded vertices. 



More precisely, define the state vector at time t to be 

s(t) = K(t) ... viit-ll) ■■■ VkH) ■ ■ ■ VKit - /^)]^. (4) 

This is the collection of all vertices at times < t, that interfere with, or are interfered with 
by vertices at times > t. Figure [3] illustrates which vertices are included in the state vector. 
As each Vj(t) takes on one of two values, the state space consists of at most 2^j=-^'^^j^^^ 
possible states. Some states may be infeasible because two of their vertices are connected 
by an edge. Thus define the state space as the space of all feasible states 

S = {Si, . . . ,S|5|} 

with each Sj G {0, 1}^j=i^'3"'"^\ Each state in S corresponds to an independent set in the 
subgraph made up of vertices in s{t). Thus there are typically far fewer than 2^j=i'''j^^'' 
states. A second example is given in figures H] and [H In this example there are a total of 
28 states as shown in figure [5l 

For notational convenience we label the elements of Sj as such 

^ _ \M M) AK,o) {K,t-i*j,).T 

Denote the set of feasible state transitions from a G iS to b G 5 by 

= {(a, b) : fe^^'^^+i) = a^^'''\ for j = 1, . . . , and A; = 0, . . . , /* - l} 
For the example of figures H] and [5] the set of feasible state transitions is 

S = {(Si, S2), (Si, Sg), (Si, Sii), (S2, S3), . . . } 

With the state space clearly defined, it is straightforward to derive the actual algorithm, 
which is akin to the Viterbi algorithm. It begins by initializing the costs to zero, that is 

cs,(0)=0, 
10 



Figure 5: An illustration of the entire state space S for the example of figure HI For each 
state, both the state vector is given and the corresponding independent set shaded. 
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for all Sj G S. Each iteration of the algorithm involves finding the minimum cost path 
entering each state. At times t = 1, . . . ,T compute 

Si*{t, Sj) = arg min Cs^(t - 1), 

Si G 5 
isi,Sj) G T 

for each Sj G S, which is the minimum cost state at time t — 1 from which we can transition 
into state at time t. Then compute 

K 

Cs,{t) = Cs^*(t,s,.)(t - 1) -^rksf''^\ 

k=l 

for each Sj G S, which is the minimum cost of a path that ends up at state s^- at time t. 
When time t = T is finally reached, compute 

s*(r) = argmincs,(r), 

Sies 

which is the optimal termination state. The optimal independent set is then found by 
working backwards. Start by setting the Tth column of vertices Vi(T), . . . , vk(T) according 
to s*(T). That is, set Vj{T) = [s*(T)](^'°) for j = 1, . . . , iT. Then set the T - 1th column 
according to Si.(T, s*(T)), that is, set Vj{T - 1) = [s*(T, s*(T))](J'°) for j = 1,...,K. 
Continue by setting Vj{T - 2) = [s*(T - 1, s*(T, s*{T)))Y^'°\ etc... 

We now briefly examine the complexity of this algorithm. As the state space consists 
of all independent sets of a subgraph defined by 2^j=i*-'j'^^^ vertices, given a set of delays 

it takes 0(2^^='^^''^'^^^) time steps to enumerate. Once this is done, the algorithm 
takes 0(T\S\) time steps to solve problem ([3]). Roughly speaking, l'-^ = 0{L), and |iS| = 
0{LK log LK). Thus the algorithm takes 

0{T LK log LK) +0{2^^) 

time steps to compute the optimal independent set. 

As the state space defined in equation H] is finite, for large T the optimal independent set 
will have a periodic form with period less than or equal to the number of states |iS|. Thus, 
if there is no restriction on how large T can be, once the period of the maximal independent 
set is found, we can simply set T equal to it, without compromising optimality. In this 
case, the entire problem can be solved in 

0{LK log 1X2^^^"^^^^) + O (2^^) = 0{{LKf^+HogLK) 

time steps. 

4.2 Bandwidth scaling 

Theorem 13.11 shows that if the bandwidth scales sufficiently quickly with K, the spectral 
efficiency per user can be made to vanish arbitrarily slowly. A natural question to ask 
is whether it is necessary for the bandwidth to scale with K, in order for this desirable 
property to hold. The following converse result establishes that this is indeed the case. 
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Theorem 4.2. // the bandwidth scales sufficiently slowly with K such that 



\ogW 
lim 

K^oo log 



K 



logK 



then 

with probability one. 

As the total spectral efficiency 

K 



Ito = (5) 



theorem 14.21 implies 



<^log2 (l + max|/^,,p^j 



log (E.=i R^ 
lim ^ — ^ = 0, 

/sT^oo log A 



which is equivalent to \im.K^oo Ri = for almost all users i G {1, . . . , K}. 

Roughly speaking the above result says that if the bandwidth scales slower than 0{K/ log K), 
then the spectral efficiency resulting from any feasible transmit pattern will vanish as 
K ^ oo. Note there is a gap between this converse result and the achievability result of 
theorem 13. 1[ Theorem 13.11 demonstrates that it is sufficient for the bandwidth to scale like 
0{{2K{K - l))^(^-i)), but theorem shows that it is necessary for the bandwidth to 
scale only as fast as 0{K / \ogK). This establishes that the slowest possible bandwidth 
scaling lies somewhere between 0{K / \ogK) and 0{{2K{K — l))-'^'^-'^"^)). It is unclear, 
which if any of these bounds is tight. 

Proof. Remove those edges that connect vertices of different time slots in the interference 
graph, {vi{t), Vj{t')) for t' ^ t. This provides an upper bound on the independence number. 
Now consider a single column V{t) = . . . ,VK{t)} of this graph in isolation. For any 

pair of vertices in V(t), there exists an edge connecting them independently with probability 
1 — (1 — l/L)"^ (probability 1/L for each of the two possible directions). Thus the graph 
consisting of vertices V(t) and the random subset of edges connecting them, is precisely 
the Erdos-Renyi graph Gk,i-{i-i/l)^- A well known result (see for example [8]) is that 

lim - 2 



logn log(l/(l-p)) 
with probability one. Hence the independence number of Qk,2/l satisfies 

a{QK,2/L) _ 2 



lim 



\ogK log(l/(l - 1/L)2) 
< L. 
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As L is directly proportional to W , if \\mK-*oo log W/ log j^^^ = then the same limit 
applies for L and 

\oga{gK,2/L) „ 
lim — = 0. 

K^oo log K 

Now as the independence number of the interference graph satisfies Oi{Q) < Ta{QK,2/L), 
equation ([5]) follows. □ 

4.3 When is the maximal independent set maximal? 

We now turn to the problem of analysis. Ideally we would like a simple characterization 
of the size of the maximal independent set in terms of the parameters of the system, 
K, T, and the normalized cross-delays It is unclear if such a characterization exists. 
Instead we present two results. The first characterizes when the independence number is 
equal to its maximum possible value, and shows how the maximal independent set can be 
found almost instantly in this event. The second, which is theorem 13. 11 demonstrates that 
surprisingly large independent sets exist on average, when both the number of users and 
the bandwidth are sufficiently high. In this section we present the former result, in the next 
section we present the latter. What we will be revealed in this section is that the problem 
of determining whether or not the independence number is equal to its maximum possible 
value, is a group theoretic one. 

We assume in this section that the direct gains are all equal so that the optimal inde- 
pendent set is equivalent to the maximal independent set. Owing to the absence of some 
edges, the boundary of the interference graph has a slightly different structure than it's 
interior. In order to circumvent this problem, we let T — oo so that these boundary effects 
are negligible. 

Definition 4.3. The independence rate of sequence of interference graphs 



is 

HGKmW,)) = lim 



- ^lim - . 



We write \R{Qk) for short. Start with the following observation. 
Lemma 4.4. For any number of users K and any channel 

This means that for large T we can only include at most half the vertices of the interference 
graph in any feasible transmit pattern. We now ask, when is the independence rate exactly 
equal to K/21 The following result succinctly answers this question for K = 2i. Define 

^ = ^31 + ^13 + ^32 + ^21 + ^12 + ^23 

h = + ^32 + ^21 

h — ^21 + ^12 

h — ^31 + ^'l3- 
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If li 7^ 0, define 7^ to be the exponent of 2 in the prime factorization of li, that is /j = 2^^/?j 
where Pi represents the rest of the prime factorization. If Zj = then define 7j = 00. 
Similarly if / 7^ 0, define 7 to be the exponent of 2 in the prime factorization /, i.e. / = 2'^ (3. 
If / = then define 7 = 00. 

Theorem 4.5. IR(^3) = 3/2 if and only 2/71 < 72, 71 < 73 and 71 < 7, in which case 
there are exactly 2S'^*^('^''2/^'^3/^'^/^) feasible transmit patterns achieving it. 

To clarify, if for example both 71 = 00 and 72 = cxd, then the above conditions are not 
satisfied and IR(^) < 3/2. This theorem provides a necessary and sufficient condition such 
that all users can transmit half the time without interfering with one another. The most 
probable way this condition can be met is if li is an odd number, and /, I2 and 1^ are all 
even numbers. Each of these events roughly occurs independently with probability 1/2, 
hence the probability all four occur simultaneously is 1/16. Thus with probability > 1/16 
there exists a feasible transmit pattern enabling all users to transmit half the time without 
interfering with one another. There are of course other ways in which our condition can be 
met, for example, if /i is even, but not a multiple of 4, and /, I2 and I3 are all multiples of 
4, however this, and all other configurations satisfying the condition of theorem 14.51 likely 
occur with probability much less than 1/16. The proof of theorem 14.51 is given in the 
appendix. 

Example 4.6. For the channel in figure [771 we have I = 2 and li = l,l2 = 2,^3 = 
2. This means 7 = 1,71 = 0,72 = 1,73 = 1, so a feasible transmit pattern achieving 
independence number 3/2 exists. As gcd(/i, ^2/2, /3/2, Z/2) = 1 there are only two feasible 
transmit patterns : the first is shown in the figure as a sequence of shaded vertices, the second 
is obtained by complementing the transmit pattern, i.e. unshading the shaded vertices, and 
shading the unshaded ones. 

How does theorem 14.51 generalize for an arbitrary number of users, Kl Define a cycle 
on the interference graph to be the indices of a tuple of edges with connecting vertices, 
that start and end on the same row. For example ((1, 2), (2, 3), (3, 1)) and ((3, 2), (2, 3)) are 
examples of cycles for K = ?>. The length of a cycle is the sum of the normalized cross-delays 
associated with it. For example, the cycle ((1, 2), (2, 3), (3, 1)) has length l'^^ + ^23 + ^si- 
Define the set of cycles containing an even number of terms as 

ye = {((ii, ^2), («2, ^3), • • • , (^2n, ^i)) : «i 7^ «2 7^ ' ' " 7^ «2n and ij G {1, . . . , K}} , 
and the set of cycles containing an odd number of terms as 

yo = {{{h, ^2), («2, ^3), • • • , (Wi, «i)) ■ h 7^ i2 7^ ■ ■ ■ 7^ Wi and ij e {1,..., K}} . 
Then we claim it can be shown that 

Claim 4.7. \R{Qk) = K/2 if and only if the exponents of 2 in the prime factorizations of 
the lengths of all cycles containing an odd number of terms, are the same, and this exponent 
is strictly less than the exponent of 2 in the prime factorization of the length of every cycle 
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containing an even number of terms. That is, the exponent of 2 in the prime factorization 
of 

is the same for all Yo G 3^oj (^iT'd this value is strictly less than the exponent of 2 in the 
prime factorization of 

for any e 

5 Achieving Non- Vanishing Spectral Efficiency 

We now prove theorem 13.11 by presenting a construction with expected spectral efficiency 
that can be made to vanish arbitrarily slowly as — oo. First, a high-level overview of the 
proof. The idea is to construct a transmit pattern that has close to 0(1) independence rate 
as i^' ^ oo. At the heart of the transmit pattern is a generalized arithmetic progression. 
If the bandwidth scales appropriately with the number of users then this progression will 
have desirable interference alignment properties, but care has to be taken in constructing 
a transmit pattern out of it. In particular, the progression will be very sparse, meaning 
that many identical versions of this progression must be interleaved, each with a different 
timing offset. The trick to making the analysis work is to use a randomization argument 
to show that a good set of offsets exists. 

Proof, (of Theorem 13. ip First some notation that will be used throughout the proof. Let 

N = K{K -1) 

and 

We use © to denote addition modulo L. Let 

^=1 a,,/,,:{a,,h^,G{0,...,Ar-l}^l. 

This is the set of all linear combinations of the cross-delays (not normalized) with integer 
coefficients ranging from to A^. Define 

A 
a=l 

See figure [6] for an illustration. Each user transmits one data symbol at each time slot in 
the set 

oo 

X ^ [j{S + kL). 

k=0 
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p rr n 



i i i 



n 



i i i 







r+ nil mod L 










r + m3 mod L 







r + m2 mod L 



time slots 



L-1 



Figure 6: The construction S is formed by interleaving a sufficient number of generalized 
arithmetic progressions T with random offsets rrij. The colored bars indicate time slots 
during which data symbols are sent. 



The above construction corresponds to concatenating data blocks S of length L. The 
modulo L addition used in constructing S and T, ensures a seamless transition at the block 
boundaries. This is illustrated in figure [71 The construction is defined in this seemingly 
convoluted way in order to make the analysis simple and elegant. However there is an easier 
way of conceptualizing this construction: take multiple copies of the generalized arithmetic 
progression |^i<j^j<^ o^ijhj '■ ^ {0, . . . , — and throw them down on the 

infinite time axis with offsets mi, . . . , rriA, L + mi, . . . ,L + rriA, 2L + mi, . . . ,2L + mA, ■ ■ ■ ■ 
Although this construction is periodic with period L, locally, the offsets of these progressions 
will appear as a Poisson process with intensity A/L = As there are A^^ points in 

each progression, the density of points in X will be and hence the spectral efficiency 
will go to zero with K \l]^e 1/{K{K - l)Y . 
We will show that there exists a choice of 



mi. 



,mA) e {Q,...,L-iy 



such that the expected spectral efficiency of this scheme approaches ([2]) as — oo. More 
specifically, we show that for the above construction, at each receiver the expected fraction 
of time slots containing a data symbol but no interference is large. Each such data symbol 
is then able to convey log2(l + \hii\^PSD /Nq) bps/Hz of information and the expected 
spectral efficiency achieved by the scheme for user i is the fraction of such time slots 
multiplied by log2(l + PSD /Nq). 

Since our construction X consists of a concatenation of identical blocks of length L, we 
analyze its performance over a single block extending from time slot to L — 1. At receiver 
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S S+L S+2L 







INN MM N 1 II 




III N 1 II 




INI III Nl II 



L time slots 2L 3L 

Figure 7: The construction X is formed by concatenating blocks of S. The modulo L 
structure of S ensures a seamless transition from block to block. 



i the set of time slots containing interference is 

K 

U {S®li,). (6) 

Denote the number of time slots at receiver that contain a data symbol from trans- 
mitter but no interference by Si = |{t G (5 © Then conditioning on the cross 
delays we have 

^S, = j^^^ (7) 

Define 

[0, otherwise. 

That is, s{k) = 1 if a transmission takes place at time slot k, and zero otherwise. Similarly 
define 

1, ifkeJ^i 
0, otherwise. 



m 



That is, fi{k) = 1 if there is interference during time slot k at receiver i and zero otherwise. 
Conditioned on the cross delays. Si is the correlation function between the set of transmit 
times S and the set of interference free times J-'f , evaluated at an offset of la, specifically 



L-l 



5',(y = 5^s(A; 



k=0 
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Thus 



L-l 



hi=0 
^ L-l L-l 

^ L-l L-l 

fc=0 lii=0 

^ L-l L-l 



lii=0 



k=0 



where in the above sequence of equations we have used the identities \S\ = J2k=o 

l-^il = J2k=o fiW- Substituting back into equation ([7j) we find the fraction of time slots at 

receiver i containing data but no interference is 



K-i) r I ^ 



IK(K-I) 



\T- 



(8) 



The above expression makes intuitive sense as if we uniformly select a time slot at random 
from {0, . . . , L — 1}, then conditioned on the the quantity \S\IL is the probability 

this time slot contains a data symbol, and 1 — \Ti\lL is the probability it does not contain 
an interference symbol. We now compute appropriate bounds on the terms \S\ and \T^. 
From equation ^ we have 



u 

A 

u 

a=l ijj^i 



\J{ma®T) 



a=l 



But 



|J(r © kj) C \ ai.kj ■ {ai,},^, e {0, . . . , iV}^ i 



19 



and this set has at most (A^ + 1)^ elements. This is the interference ahgnment property. 
Hence 



^^4 



< 



< 



[j{r © 

A{N + 1)^ 
L 

J\fN+e 



= AT-M 1 + _ 
V N 

< eN-\ 



N 



We now bound |«S|. We first show that |T| = almost surely as — > oo. In order to 
have |T| < A^^, there must exist two sets of coefficients {oLij}i^j ^ {c^ijli^j both elements 
of {0, . . . , AT - 1}^, satisfying 

This is equivalent to requiring there to exist some {aij}i^j G {— A^ + 1, . . . , A^ — 1}^\0 
satisfying 

Oiijhj = 0. 

Using the union bound we have 

Pr (|r| < AT^) = Pr ( 3{aij}i^j e {-A^ + 1, . . . , A^ - 1}^\0 s.t. ^aijkj = | 

{a^,h^,€{-N+l,...,N-l}'^ / 

As conditioned on all cross delays other than /12, there is at most one value of /12 that 
satisfies ^j^j c^ijlij = 0, and the cross delays are uniformly distributed over {0, . . . , L — 1}, 
we have 



Pr (|r| < AT^) < 

{aij}ij,j&{-N+l,...,N-l}N 



1 

L 



< 



L 
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overall density of points = 0(2"'^) 



N distinct points 




0(N'^) time slots 



0([2N]^) time slots 



Figure 8: An illustration of the generalized arithmetic progression T without the modulo 
L wrap around. There are points spread over 0{{2N)^) time slots, however almost all 
points are concentrated at the center in a width of 0{N^). Hence the density of points is 
0(2- 



as N oo (or equivalently, as K ^ oo). 

At this point it should start to become clear why it is that in this particular construction 
the bandwidth must scale like {2N)^ . From the above calculation we see that in order to 
make all the points in the generalized arithmetic progression T distinct, we require the lij to 
be large. How large? It may seem that as there are A^^ integers in T, the minimum being 
and the maximum being roughly the same order as lij, we require 1^ = 0{N'^). However, 
the structure of the generalized arithmetic progression is such that the bulk of its points 
are concentrated around the center, such that we actually require at least 1^ = 0((2A^)^) 
to separate these center points out, as the above calculation shows. But such a large 
order of kj makes T very sparse, in fact if Uj = 0{{2N)^) then T's density is a mere 
0{2~^). So to fill in the gaps we interleave multiple sequences T. How many? 0(2^). In 
general if kj = 0{L) then we must interleave 0{L/N^) = A sequences. This explanation 
is illustrated in figure M 

We now use a probabilistic argument to demonstrate the existence of a good choice of 
(mi, . . . , rriA)- Let nia ~ i.i.d. f/({0, . . . , L — 1}). Let 

t^k) A / 1' if ^ e r 

[0, otherwise, 
otherwise. 



+ ■■■ + tAik)il - tA-l{k)) X ■ ■ ■ X (1 - U{k)). 

This expression says that k E S ii k & mi © T, or if ^ mi © T but k G m2 © T, or if 
k ^ mi © T and k ^ m2 © T but k G m-^ © T, etc... We can write the above expression 



and 

Uk) 4 { J; 

Thus taik) = tik © ma). Then we can write 
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alternatively as 

s{k) = h{k) + (1 - hikmt2{k) + (1 - hikmtsik) + ...)) 

= t{k © mi) + (1 - t{k © mi)){t{k © ms) + (1 - t{k © m2)){t{k ®mz) + ...)) 
Then taking the expectation over the distribution of mi, ... , m^ 



E|5| 



1 



A:=0 
L-1 



= yE V © mi) + (1 - t{k © mi)){t{k © ma) 



fc=0 

+ (1 - t{k © m2)){t{k © mg) + . . . )) 

L-l L-1 L-1 



X] ■■■^^t{k®mi) + {l-t{k® mi)){t{k © ma) 

+ (1 - t{k © m2)){t{k © ma) + . . . )) 



m^=0 mi=0 fe=0 



^ L-1 L-1 L-1 

= TaTT 5Z 5^ ■ ■ ■ 5^ © mi) + (1 - © mi)){t{k © ma) 



fc=0 mA=0 mi=0 



+ (1 - t{k © m2)){t{k © ma) + . . . )) 

L-1 L-1 L-1 

4TrEE---El^l + (^-|^IM^®^2) 



fe=0 mA=Q m2=0 



+ (L - |r|)(l - t{k © m2)){t{k © mg) + . . . )) 

^ L-1 L-1 L-1 



fe=0 7nA=Q m2=0 



l/j-l ^ — 1 / \ a 



+ (L-|T|)2(t(fc© ma) + ...)) 



El 



a=0 



1 - 1 - — a.s. 



1 _ g Z, 
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Substituting back into equation ([8]) 




L 



(1 - eA^-^) 



as N oo (or equivalently K oo). As each data symbol that is received without inter- 
ference is capable of reliably communicating log2(l + I/im^PSD/No) bps/Hz, the expected 
spectral efficiency of each user i goes to 



The proof of theorem 13.21 is a straightforward extension of the previous. 

Proof, (of Theorem 13. 2p Assume without loss of generality that argmax^^ Iha^dl"^ = 1 for all 
users i. If each receiver i treats physical paths 2,3, ... ,D from transmitter i as interference, 
then the received signals in the K user D path interference channel are statistically identical 
to those of the DK-user LOS interference channel. Thus the achievability result of theorem 
13.11 carries over to the D path channel with K replaced by DK, and replaced by 

maxrfg{i^...^£)} \hii^d\'^. □ 

6 Frequency Domain Interpretation 

In this section we reconcile the time domain version of interference alignment presented in 
this paper, and the frequency domain results of [5j. Specifically we show how the three-user 
construction of |Sj has a simple time domain structure for the LOS interference channel. 

To begin, we need to transform the LOS model into the frequency domain. For this we 
use an OFDM architecture summarized in figure [H Transmitter j has a stream of complex 
data symbols to send {xj[0], . . . } to receiver i. These are broken up into blocks of 
length n. Consider a single block denoted x = [xj [0], . . . , Xj [n — 1]]"^. To send this block the 
transmitter computes the M- length vector xj = V^Xj, where Vj G C^^^" is an encoding 
matrix to be specified later. Let /max = maxj^j kj. To send Xj, tx j computes its IDFT and 
appends a cyclic prefix of length Zmax- Each receiver removes the cyclic prefix and computes 
the DFT. Specifically 



1 



-log2 (l + |/i.i|2pSD/No) 



{K{K-1)) 



as K 



oo. 



□ 




and 
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Figure 9: Illustration of the OFDM architecture used to reconcile the time and frequency 
domain versions of interference alignment. 



where Fmxm is the M x M DFT matrix, 



/ 1 

1 



MxM 



1 

gi27r/M 



J2-K-2/M 



\ 



I gi27r-2/M gi47r-(n-l)e' 
I gi27r-(M-l)/M gi27r-2(M-l)0 



1 \ 

gi27r-(n-l)/M 

gj27r-2(n-l)/M 
^i2ir-(M-l)(n-l)/M 



The result is a length M sequence 

K 



-j2Tr{fcnj+klij/M) 



for A; = 0, . . . , M — 1. Let Oij = kj/M. Then in matrix form 



K 



where the link matrices are 



jr.. — h-.p-j'^-^fcnj 



( 1 



g-i27r-eij 



g-i27r-2eij 



V 



g-j27r.(M-l)0i, 



(9) 



Choose W sufficiently large such that the lij are distinct. Note that M needs to be much 
larger than /max in order for the overhead from the cyclic prefix to be small. Also, M 
must not have any of the Uj as divisors, else the channel matrices will lose rank. Let 

H = /ir.ie^'2'^/=^».H and T = Hi2H2/H23H32^H3iHi3\ Let w = [1 ••• 1]^. Choose the 
encoding matrices as follows 

V2 = TV 

V3 = H]^3 H12V 
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where 
and let 



V = [w Tw TV ■ ■ ■ T"- V 



n = (M- l)/2. 

where M will be chosen to be an odd number. This is the three user construction of [5], 
but note the channel matrices Hj^ do not consist of M independently faded tones. Rather, 
all tones are derived from a single parameter lij. Define 

= 612— 621 + 6*23 — ^'32 + ^'31 — 6'i3 

1 = h2 ~ hi + ^23 ~ ^32 + ^31 ~ ^13 



Then 



V 



/ 1 

1 
1 



1 

gi27r-20 



V 



I gj27r-(Af-l)e gj27r-2(M-l)0 



gi27r-(n-l)e 
3j27r-2(n-l)0 



J2n-iM-l)(n-l) 



(10) 



Lemma 6.1. // M is prime the columns of V are a permuted subset of the columns of 

V = Fmxmt^im (11) 

where tti^m is an M x n permutation matrix, i. e. each column of tt^m is a unique column 
of ImxM- 

Proof Consider the matrix element V(2, k) = e^"^^^^/^ = e^^^'^*^' ""'^ ^)/^^/ for 
some A; G {0, . . . , M — 1}. Let M be prime. Then the set 

gi27r/M^ gi27r-2/M^ _ _ _ ^ gi27r-(M-l)/M| 

together with the multiplication operation forms a group. Thus each V(2, k) corresponds 
to a unique Fmxm(2, /c') for some k' G {0, ...,M — 1}. Now observe that V(j, A;) = 
V(2,fc)(^-i) and FMxMU^k') = Fmxm{2, k'Y^''\ Thus FMxM{j,k') = Yuxuij^k). In 
other words each column of V corresponds to a unique column oIFmxMi which establishes 
the result. □ 

Lemma 16.11 enables us to write the encoding matrices V, in a revealing form. Define 



Fi — Ho]^ H32T 



To ^T 



1^3 — H12. 



(12) 
(13) 
(14) 



Then 



^^max X {iV/ /max) "^^max X ^ma 



F M X Af r j F M X Af 71" « , M X J . 
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Examining the above expression reveals that the encoding operation for tx j corresponds 
to transmitting consecutive data symbols / time slots apart, but cyclicly wrapped around 
such that roughly half of all time slots contain data symbols and no two data symbols share 
the same time slot. As the operation Fl_,j^^jrjFMxM corresponds to delaying the input 
stream, the entire transmission sequence is just offset by this amount. 

Based on equations f|T2|) - f|T^ we can define the delay dj for user j's transmission se- 
quence as 

di = 2/12 — 2Z21 + 2/23 ~ ^32 + ^31 ~" 2Zi3 
d2 — ^12 ~ hi + ^23 ~ ^32 + ^31 ~ ^13 
ds — /12 ~ ^13- 

Then we see that transmitter j will send its first symbol Xj [1] in time slot dj mod M, its 
second symbol xj [2] in time slot dj + / mod M, its third in time slot dj + 21 mod M, etc... 
The last symbol will be sent at time dj + {n — 1)1 mod M. More precisely transmitter j 
sends 

~ r 1 / if m = kl + dj mod M for some k 

"'^ \ 0, otherwise 

at times m = 0, 1, . . . , M — 1. 

The manner by which this transmission scheme achieves interference alignment is now 
simple to understand. Tx 1, transmits it's first data symbol a;i[l], such that it arrives at 
rx 3 at the same time as tx 2's first data symbol X2[l]. Tx 2 transmits a;2[l] such that it 
arrives at the same time as xsfl] at rx 1. Tx 3 transmits a;3[l] such that it arrives at rx 2 
at the same time as Xi[2], etc... See example 16.31 and figure [TOl 

From figure [TU] it is clear how decoding should be performed — receivers merely decode 
each data symbol by looking at the time slot in which it was received. But let us reconcile 
this with the decoding methodology of where the received sequence is passed through 
a ZF equalizer. This corresponds to projecting the vector onto the subspace orthogonal 
to the interference. At the first receiver 

yi = ( Hn Vi H12U ) ^1 ^ + zi 

where 

U= [wTwT^w ■■■ T"w]. 

and Xg represents a combination of interfering symbols from the 2nd and 3rd users. It is 
straightforward to see that we can write the space orthogonal to U as 

U'^ = Fmxmt^Im 

where '7ifj^,j is a permutation matrix orthogonal to tti^m in the sense that t^i = 0. 

Simply put, U'^ is a matrix whose columns are those columns of Fmxm that are not present 
in V. Thus the first receiver computes xi = G*yi where 

Gl = (U'=*Hi2^HiiVi)-iU^*Hi2^ 
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If we write out the decoder in detail 

(^?,M*FMxA/Hi2"^HiiriFMxAf7r/,A/) V;';a,j*FX^xmHi2^Fmx A/ ( OMximax ^A4xM ) Yj 

we see that the first few operations correspond to removing the cychc prefix, delaying the 
resulting stream by I12 and then selecting the interference free subset of this. The last 
operation is in general undefined, as the matrix we invert may not be full rank. However 
in certain scenarios the matrix equals identity and is then invertible. 

The reason for this phenomenon is that the last operation corresponds to recovering 
the data from the interference free subspace, which may contain fewer than n dimensions 
of Xj . In the case of many independently faded OFDM sub-channels such as is assumed in 
[5], one could always project the subspace spanned by the data onto the interference free 
subspace without losing information, however for single path channels each data dimension 
is either orthogonal or overlapping with an interference dimension. 

The scenario in which the above decoder is well defined (i.e. data and interference 
subspaces are orthogonal) is given by the following condition. 

Lemma 6.2. // 

^11 = di + nl mod M 

^22 = d2 + nl mod M 

hs = + nl mod M 

then the interference subspace is orthogonal to the data subspace. 

Thus if each of the direct delays takes on a single, specific value, the signal space will be 
orthogonal to the interference space at each of the receivers and we will be able to decode 
all data symbols. 

Example 6.3. Suppose the link delays are In = 0, I22 = 5, I33 = 9, I12 = 5, /13 = 4, /21 = 4, 
= 6, /31 = 5, I32 = 4. Then rfi = 7, c?2 = 4 and d^ = 1. Also I = 5-4 + 6-4 + 5-4 = 4. 
Note that for illustrative purposes the direct delays have been precisely chosen such that 
the data and interference subspaces are orthogonal. Choose the data block length M to be 
the prime 13, and use a cyclic prefix of length 9. The total length of cyclic prefix plus 
data block is 21. Then tx 1 will transmit its 0th data symbol, namely Xi[0], in time slot 
di mod M = 7 mod 13 = 7. It's second data symbol xi[l] will be transmitted in time slot 
I + di mod M = 4 + 7 mod 13 = 11. It's third data symbol xi[2] will be transmitted in 
time slot 21 + di mod M = 15 mod 13 = 2, etc.. These data symbols will arrive at rx 2 
delayed by /21 = 4 time slots. Thus xi[0] will appear as interference at rx 2 during time slot 
7 + 4 = 11, Xi[l] will appear as interference during time slot 11 + 4 = 15, etc... Similarly 
tx 1 's data symbols will arrive at rx 3 delayed by /31 = 5 time slots. Similarly one can do 
the same computation for tx 2's and tx 3's data symbols. 

The details are given in figure {Wi In part (a) of the figure the cyclic prefix has been 
omitted for illustrative purposes. It is incorporated into the picture in part (b). The red 
shaded boxes contain data symbols, the grey shaded boxes contain interference symbols. 
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Figure 10: Illustration of 3 user interference alignment scheme of [5] in the time domain. 
See example 16.31 for a description, (a) Received sequences with the cyclic prefix omitted, 
(b) Received sequences incorporating the cyclic prefix. 
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Notice the interference alignment property manifests itself as an overlapping of interfering 
data symbols. The shaded, but unlabeled boxes represent symbols belonging to the next 
OFDM block. The black box outlines those time slots that are used for decoding. The 9 
time slots prior to these are discarded when the cyclic prefix is removed. The 9 unboxed 
time slots to the right of the black outline will also be discarded, but during the next OFDM 
block. 



6.1 K-user channels 

In the previous section we demonstrated that for three-user LOS channels, the frequency 
domain scheme of [5] has a simple analog in the time domain that works well when the 
block length is chosen to be a prime number, and the direct delays take on particular 
values. One would imagine that an analogous scheme for K > ?> users would therefore 
also exist and work well. This is not the case. In fact for the LOS interference channel 
with more than three users, the alignment scheme of [5] has various shortcomings which 
result in it achieving zero degrees of freedom in total. Our construction (in section [5]) 
is inspired by the use of a generalized arithmetic progression in |.5j, but circumvents the 
schemes shortcomings by: 

1. Truncating the generalized arithmetic progression appropriately. 

2. Interleaving many replicas of the progression, with random offsets. 

3. Scaling the bandwidth with K. 

For example, to understand why truncation is necessary, recall the precoding matrices 
used in for the i^-user channel in [3] are of the form 



S,B 



for j 



2,. 



B 



K, where the matrix B is composed of the column vectors in the set 
/ 



\i,je{2,3,...,K},i^j,{i,j)^{2,3) 



w : aij e {0,1, ... ,n - 1} 



and Sj = H13H23 H21. Observing equation Q), write the link matrices in the form 



H, 



where 



/ 1 



-j2Tv/M 



-j2-K-2/M 



'j2-K-{M-l)/M 
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Then we have 

B = |z^»j6{2.3.....K},i7^i.{«.i)7^{2,3)°>.'»jw : \/aij G {0, 1, . . . ,n - 1}| . 

where kj = —kj + 1^ — lij + lis ~ hs + hi- For sufficiently large n we will be able to find 
many pairs ({ajj}, {o^ij}) such that 

Thus the precoding matrices will loose rank as many of their columns will be repeats of 
previous ones. This phenomenon of repeated elements is common in generalized arithmetic 
progressions over integer fields. See for example [9]. How much rank will be lost? Observe 
that the largest exponent of Z in i3 will be no greater than {n~l){{K —1){K — 2) — 1) maxlij. 
But there are n^^'^^^^'"^^"^ columns in B. Thus as n — * oo the rank of B will scale only 
like 0{n) due to repeated columns, whilst the dimension of the space scales like roughly 
0(n(^-i)(^-2)-i). Hence the total de grees of freedom goes to zero unless the progression 
is truncated. 



6.2 Bandwidth Scaling Revisited 

The final issue we address in terms of reconciling time and frequency domain interpreta- 
tions, is that of bandwidth scaling. We now demonstrate that bandwidth scaling is required 
in the scheme of [5J when the physical channel model is brought into the picture. 

Theorem 6.4. In a multipath fading channel with L taps, if the bandwidth satisfies 

y logW^ 

K^o. \og{{K -l){K-2)- l){K-l){K-2)-Z ^' 

then the total degrees of freedom achieved by the K-user interference alignment scheme of 
1^ goes to zero as K ^ oo. 

This means that the bandwidth must scale at least as fast as 0{{{K — 1){K — 2) — 
l)('^-i){^-2)-3^ ^j^ich is roughly the same scaling that is required in theorem 13. namely 

Proof. For a general multipath channel with L taps, the link matrices are of the form 

L-l 



Using the commutativity of the diagonal Hjj matrices we can write 
^3 = Si I W HiiHijHss 

Vije{2,3,...,/^},j7^i,(j,i)7^(2,3) 
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where the matrix C is composed of the column vectors in the set 



C = <j I n (H,iHi,H23)"-"-(H,,Hi3H2i)"- I w : 

^ijG{2,3,...,K},i^j,ii,j)^i2,3) 

aij e {0, 1, . . . ,n - 1} 

In [5] the minimum scahng of n with K required is 

hm > 0. 

K^oolog{K -1){K -2) - 1 

Each of the Hij matrices is a polynomial of degree at most L — 1 in the matrix Z. Thus 
each column of C is a polynomial of degree at most 6n{L — 1){{K — 1){K — 2) — 1) in the 
matrix Z. This means the maximum rank of C is 6n{L — 1){{K — 1){K — 2) — 1) + 1, as 
any c polynomials of degree < d that are in general position, are linearly dependent for 
c > d+1. The total number of rows in C however, is at least n^^~^^^^~'^^~^. Thus the total 
degrees of freedom is no more than 

Qn{L - 1){{K - 1){K - 2) - 1) 



^{K-l){K-2)-l 

which goes to zero as ^ oo, unless L (and hence W) scales like 

K^oo \og{{K - 1){K - 2) - l){^-l)(^-2)-3 ^ 



□ 



7 Discussion and Conclusion 

We demonstrated in section 14.21 that if the bandwidth scales sub-linearly in then the 
independence rate of the interference graph goes to zero as K ^ oo, and in section [5l 
that if the bandwidth scales like 0{{2K{K — l))^*^^"^)), then the independence rate scales 
arbitrarily close to 0{K). This brings us to an interesting open question. How does the 
independence rate scale in the intermediate regime where 0{K/ logi^') < W < 0{{2K{K — 
What about the D-path channel? 
There are several variations of the LOS interference channel for which the interference 
graph techniques discussed in this work are applicable, and for which further study is 
warranted. These were alluded to earlier. For instance, the partial connected interference 
channel is a more accurate model of an extended wireless adhoc network. It is possible that 
for such channels, a bandwidth scaling much less than 0{{2K{K — l))-^*^-^"^^) is sufficient 
in order for the independence rate to scale arbitrarily close 0{K). It is not clear how one 
would approach the problem of showing this if it were true, or disproving it otherwise. An 
interference channel with one dominant path per link and several sub-dominant ones, is also 
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an interesting candidate for investigation. This D-path channel is commonly encountered 
in practice. It is possible that an optimization problem similar to [31 but allowing for signal 
and interference to overlap, can be formulated for this scenario. It would be interesting 
to study whether time-domain based interference alignment can provide gains in this case. 
Presumably, for fixed W, as the dominance of one physical path over the others diminishes, 
so too will the performance improvement. 

Lastly we discuss some interesting fringe benefits associated with the communication 
schemes presented in this work. Whereas interference alignment in the frequency domain 
requires coding over very long blocks, which results in substantial delay due to the necessity 
of buffering data symbols at the encoder and received symbols at the decoder, no such delay 
is required for the time-indexed interference graph techniques detailed above. Data symbols 
are transmitted as soon as an appropriate time slot is reached, and detected when received. 
The only delay incurred is that stemming from the use of an error correction code. In the 
same respect, the encoding and decoding complexity are greatly reduced. Thus delay and 
complexity issues are non-existent. 

Subspace conditioning issues are also non-existent. Interference alignment in the fre- 
quency domain, although performing well at very high-PSD, suffers at moderate PSD if 
the data and interference subspaces are close to one another. Time domain interference 
alignment techniques are free from this problem as the data and interferences subspaces 
are orthogonal by design. 
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9 Appendix 

9.1 Proof of Theorem 14.51 

First suppose / 7^ 0. From the interference graph form I infinite chain graphs Q'q, . . . , 
These graphs will be functions of K and but for notational brevity we omit this 

notation. The ith chain graph has vertex and edges sets 

K = K(o),t':(i),...} 

^; = {e:(o),e:(i),...}. 

This is an undirected graph with edge e[{j) joining vertices fj'(j) and f^(j + 1). We now 
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(a) interference graph (b) one cycle of interference 
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V3(0), ... ,V3(7) 




Vu(0) v„(2) 
& v,(6) 



(c) chain graphs 
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Figure 11: (a) A segment of the interference graph. Each row of vertices corresponds to 
the transmission opportunities for each of the users. The shaded vertices correspond to a 
feasible transmit pattern achieving an independence rate of 3/2. A few of the vertices are 
labeled with their equivalent vertices in the chain graphs. In this example the normalized 
cross delays are /21 = 0, I'^i = 1, = 2, l'^2 = 0; ^'13 = 1 ^23 = ~2. Thus / = 
+ 1 + 2 + + 1 — 2 = 2 meaning that each cycle of interference moves two time slots 
to the right as illustrated in (b), which shows a single cycle from the interference graph, 
containing all six directed edges, (c) The corresponding chain graphs. As / = 2, there 
are two chains. A pair of twin vertices connected by a dashed line in the chain graphs, 
correspond to the same vertex in the original interference graph. 



associate the chain graphs with the interference graph. Let 



v-{6k) =v3{kl + i) 
v'-{Qk + l) = vi{kl + i + 
v'i{6k + 2) = v3{kl + i + 
v'i{6k + 3) = V2{kl + i + 
v-{Qk + A) = vi{kl + i + 
v'i{6k + 5) = v2{kl + i + 



hi, 
II 

'■31 
il 

'31 
II 

'31 
ii 

'31 



+ /^g + 1^2. 
+ /;^3 + /gg 

+ /;^3 + 



+4 
+ /'+/; 



12) 



where the Vk{t) are the vertices in the original interference graph. See figure [TT] for an 
illustration. Note that the mapping from the interference graph to the / chain graphs is 
not one-one. In particular, each vertex in the interference graph is associated with a pair 
of vertices in the set of chain graphs. Paired vertices are called twins as they correspond to 
the same vertex from the original interference graph. We now think of a feasible transmit 
pattern as a collection of vertices from the set of chain graphs. However, note that if a 
feasible transmit pattern results in a particular vertex (from the set of chain graphs) being 
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included in the independent set, its twin will also be included. Likewise if a feasible transmit 
pattern results in a particular vertex being excluded, its twin will also be excluded. The 
key to characterizing the set of channels for which an interference rate of 3/2 is achievable 
lies in understanding which pairings are favorable, and which are not. 

The pairings can be succinctly described by the following three equations 

vlW = vii^^^^Mk + [h/l\) + 2) (15) 
vliQk + 1) = vli^^^Mk + [h/l\) + 4) (16) 
v',{6k + 3) = vli^^^Mk + [h/l]) + 5) (17) 

for 2 = 0, — 1 and A; = 0, 1, ... . 

In order to achieve an independence rate of 3/2, half of all vertices must be included 
in the transmit pattern. Denote the transmit pattern by T. Because in each chain, all 
neighboring vertices are connected by an edge, this is only possible if in each chain, every 
second node is included in the transmit pattern. For each chain there are two ways of doing 
this, either v[{2k) G T for all k, or v[{2k + 1) G T for all k. Let q denote the phase of the 
ith chain. If the former condition holds, we say the chain is in phase and write q = I. If 
the latter holds we say the chain is out of phase and write q = O. In the entire graph there 
are only 2' combinations we need to examine, corresponding to all possible inphase/out of 
phase assignments for the / chains. A feasible transmit pattern achieving independence rate 
3/2 exists if and only if each chain admits an I or O assignment and the assignment of I's 
and O's to the / chains does not violate conditions (IT^ - (IT7|) . Thus we wish to characterize 
those channels for which such an I/O assignment can be found. 

At this point we consider an example. Suppose / divides li. We claim an independence 
rate of 3/2 is not achievable. To see this argue by contradiction. Assume that cq = /. As / 
divides Zi, we have li mod / = and condition f|T6|) pairs vertex vq{1) with vertex fo(6/c' + 4) 
for some integer k' . But these vertices lie an odd distance apart on the same chain, so 
working backwards we see that we must have Cq = O, a contraction. So suppose instead 
that Co = 0. Using the same logic as before we arrive at cq = /, again a contradiction. 
Thus the independence rate is less than 3/2. 

From this example we see that condition (fTSjl tells us if cq = / then we must also 
have C-igmodi = h as vertices v'q{^) and + 2) are an even distance apart. 

Continuing this logic we see that we must also have c.^^mod/ = O, and c^i^^aodi = I- We 
can also conclude something else, as C-l-^ modi = O, we must have c_2«i mod/ = I condition 
f fT6|) . Continuing further, we must have c_2ii-/2modz = I condition f|T71) and so on. 

By this point it should be clear that conditions (|T5|) - (|T7j) are satisfied if and only if for 
all integers ^1,^2,^35 

C2fci/i+fc2/2+fc3'3mod/ 7^ C(2k^j^l)h+k2l2+kzhmodl (18) 

Let Vil) denote the group consisting of integers {0, . . . , / — 1} together with the addition 
modulo / operation. Consider the set of chains C2A:i/i+fc2/2+A:3/3mod/ for all integers /ci, /c2, k^. 
This set forms a subgroup of V{1) with generator gcd(2Zi, Z2, ^3)- We denote this subgroup 
by 'Pgcd{2/i,/2,/3)(0- It has gcd(2/i, ^2, ^3, — 1 cosets other than itself. The set of chains 
C{2ki+i)h+k2i2+k3hmodi for all integers ^1,^2,^3, forms coset number /i modgcd(2/i, ^2, ^3, 0- 
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But as 

2/imodgcd(2/i,/2,/3,0 = 0, 

either /i modgcd(2/i, ^2, ^3, equals gcd(Zi, /2/2, ^3/2, //2), or 0. If it is zero, condition f|T8l) 
above is violated. This occurs if and only if li is a multiple of gcd(2Zi, I2, 13, 1). Alternatively 
if it equals gcd(Zi, /2/2, Z3/2, 1/2) then we can choose the phases of half the cosets, namely 
cosets 

0,l,...,gcd(/i,/2/2,/3/2,//2)-l 
arbitrarily, and still satisfy (ITS!) . For this reason we refer to the chains 

Co, Ci, Cgcd(/i,/2/2,Z3/2,V2)-l 

as seed chains. This means that there are 2^'^'^^^^''2/^''3/^''/^-' possible solutions that achieve 
independence rate 3/2. So what does it mean for li to not be a multiple of gcd(2/i, I2, 13, /)? 
It means that 

gcd(/i,2/i,/2,/3,0 ^gcd(2/i,/2,/3,0- 

In other words 

gcd(/i, k, h, I) 7^ gcd(2/i, I2, h, I). 

It is shown in lemma I^TT] that the above inequality is equivalent to having 71 < 72, 71 < 73 
and 7i < 7. This establishes theorem 14.51 for / 7^ 0. 

Now suppose 1 = 0. This proof is a slight modification of the previous. From the 
interference graph from an infinite number of cycle graphs Qq,Q[, . . . . The ith cycle graph 
Q[ has vertex and edge sets 

V: = K(0),t;:(l),t;:(2),t;:(3),t;:(4),.;:(5)} 
= {e:(0),e:(l), 6^(2), e:(3), 6^(4), e:(5)}, 

where edge e[{j mod 6) joins vertices Vi{j mod 6) and f j(j+l mod 6), for j = 0, 1, 2, 3, 4, 5 
Notice that for / 7^ we created a finite number of chain graphs, each with an infinite num- 
ber of vertices, whereas for / = we create an infinite number of cycle graphs, each with a 
finite number of vertices. We now associate the cycle graphs with the interference graph. 
Let 

^K4) =t^l(^ + /31 + ^i3 + 4 + 4) 

v[{h) = V2{i + I'^i + I'n + 4 + ^21 + I'u) 

where Vk{t) are the vertices in the original interference graph. As before, the mapping from 
interference graph to the indefinite number of cycle graphs is not one-one. Each vertex in 
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the interference graph is associated with a pair of vertices in the set of cycle graphs. The 
pairings are described by the following three equations 

%'(0)=^U(2) 

We want to assign each cycle graph i a phase, either q = J or Cj = O and we need to 
find necessary and sufficient conditions for a feasible assignment. Similarly to equation [T8| 
the condition we are after is 

C2kih+k2l2 + k3l3 7^ C(^2ki + l)h+k2l2+k3l3 

for all integers fci, A;2, ^3. By this point it should be clear, based on the proof for the / 7^ 
case, that the above condition is equivalent to 

gcd(/i,/2,^3) ^ gcd(2/i,/2,/3)- 

This inequality is equivalent to having 71 < 72 and 71 < 73. Thus we see that the conditions 
for / = case are the same as the / 7^ case if we set 7 = 00. This establishes the result 
in general. 

Lemma 9.1. gcd(/i, I2, 13, 1) ^ gcd(2/i, I2, h, I) if and only 2/71 < 72, 71 < 73 and 71 < 7. 
Proof. Suppose 

gcd(/i, k, h, I) ^ gcd(2/i, I2, h, I) 
^ gcd(2^^A, 2-^^/52, 2^^/53, 2^P) ^ gcd(2^^+Vi, ^"^2, T'^Ps, 2^/5) 
^ gcd(2^S 2^% 2^^ 2^) gcd(/?i, /32, /?3, (3) ^ gcd(2^i+\ 2^% 2^^^ 2^) gcd(/?i, /?2, Ps, P) 
=^ min(7i, 72, 73, 7) 7^ min(7i + 1, 72, 73, 7). 

This can only hold if 71 < 72, 7i < 73 and 71 < 7. The proof in the opposite direction is 
identical. □ 
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